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1. INTRODUCTION 


In the exospheric description of “'>lar-wind plasma (Eviatar 
and Schulz, ^968; Perkins and Spight, 1970), individual elec- 
trons and protons execute trajectories that conserve both total 
energy (kinetic plus potential) and first adiabatic invariant 
(magnetic moment). The motion of these part'cles along ^ (the 
interplanetary magnetic field) is profoundly affected by an 
electrostatic field that arises in order to assure charge quasi- 
ne trail ty of the plasma, i .e. . to equalize the escaping fluxes 
of protons and electrons in the presence of the Sun's gravi- 
tational field (Jockers, 19^3; 1970). 

Calculation of the corresponding electrostatic potential 
cp, at a point where the electron flux is locally isotropic and 
Maxwellian except for the escape trajectories, has been described 
by Schulz and Eviatar (1972), The local potential op is there- 
by related to the pro on bulk speed u along ^ by means of expo- 
nential and error functions of -qg^p/nTg, where h is Boltzmann's 
constant, T Is the electron temperature, and q (<0) is the 
electronic charge. Observations by Feldman et a 1 , (197l»a), 
however, suggest that the solar-wind electron distribution at 
the Earth's orbit is more nearly a superposition of two quas?- 
Maxwellian distributions, related in temperature such that T 2 ~ 6T 
and in density such that N2~0.06Np Moreover, the major elec- 
tron compont. t seems to have a velocity distribution that is 
symmetric about V||»«u secx - c^ fwhere X is the local spiral angle 
cosx * r • B, and c. Is the Alfven speed) rather than about v,| = 0 
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(as exospheric theory would pre'^{ct). 

These new observational findings require that the exospheric 
theory be updated, but not necessarily discarded, in view of 
the plasrr.a-kinet ic instabilities (Forslund, 1970) to which an 
electron exosphere Is subject (Schulz and Eviatar, 1972; Perkins, 
1973), it was not really expected that the observed distribution 
would be symmetric about v^ * 0 (Hollweg, 1974). Moreover, 
much can still be learned about the consequences of such plasma 
instabilities by treating the minor (hotter) electron compo- 
nent exospherical ly. By increasing the outward bulk velocity 
of the major (cooler) electron component, for example, the un- 
stable interplanetary-plasma waves may increase the loca^ electro- 
static potential op required to enforce charge quas i -neutral i ty. 
Such a change in cp would profoundly, although indirectly, af- 
fect the heat-flow moment of the entire electron distribution. 

It is the purpose of this work to examine such plasma ef- 
fects in detail. First, by treating both electron components 
exospher ical 1 y , we hope to estimate cp locally for the hypothetical 
(but unstable) electron exosphere. Second, by treating only 
the minor (hotter) component exospher ical 1 y, we hope to esti- 
mate the increased qp that is empirically compatible with the 
effects of the inevitable plasma turbulence. The analytical 
evaluation of heat-flow moments in both cases will serve to 
distinguish the two cases in a dynamically interesting way. 
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2. TWO- COMPONENT ELECTRON EXOSPHERE 


The first step In constructing a properly normalized exo- 
spheric distribution function f(vi|,Vj^;r) is to select a closed 
surface r > rg (not necesiirily a sphere, since rg could be 
a function of solar latitude) on wh I ch cp « tpg and beyond which 
collisions are considered negligible. Since gravity Is negligible 
for electrons, the function f(vn,Vj^;r) must be even with re- 
spect to V(i for electrons satisfying the condition 

+ 2qg«pQ <0 ( 1 ) 

2 

and must vanish for vu < 0 for values of v not satisfying (1). 

An acceptable choice ( cf . Schulz and Evlatar, 1972) turns out 
to be 

2 

f(vii.v^;rg) - Y fj(v„,v^;rg) 

j = l 

- 2[e(v„ ) + e( - vn )3 ( - m^v^ - 2q^epQ)] 

2 

X ^ N j ( m^l 2tthT . ) ^exp ( - m^v^/ 2 kT . ) 
x[l + erf( ( - q^(Pg/KT.)’^2) 

- 2( - dgCPo/tTHTj /^'^^expIqgflpgATj )]■ \ (2) 


- 7 - 


whi^re 9(x) i$ the unit step function [e(x) - 1 for x»0 and 
9(x) « 0 for x<0]. if Np is the proton density and alpha par> 
tides are neglected, quasi -neutrali ty thus requires that ■ 
N 2 and that 


MpU sec % 


M I ''i ‘*''1 


• 0 


(3) 


The integral in (3) is most easily evaluated in spherical co- 
ordinateSt and the result is that 


(2TTm^u^/xT| l^^^secx * 

2 

j-1 

X [I +erf ( ( 

■ 2 i ^^^exp(q^(PQ/»Tj )]"^ (4) 


This equation determiner (- q^cpQ/*^T^) as a function of 
(2Ttm^u^/HT^ ) ^^^sec X for fixed values of Nj/Np and T 2 /Tp Rep- 
resentative contours are plotted in Figure I. The contour N 2 /Np 0 
corresponds to the case treated by Schulz and Eviatar (1972), 
who also made a minor normalization error that Is corrected 
here. 

It is instructive to define partial heat-flow moments 

^ » 2nJ J (mg/2 ) (v^ - u sec y )^f. Vj^ dv^ dvj (5a) 

- eP 0 
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and 


q(J) 


2tTf J (m^/2)vf(v| - usecx)^jVx (5b) 

T.« 0 


for the “conduction of thermal energy" along For Maxwellian 
distributions truncated as In (2), the moment integrals are 
most easily evaluated in spherical velocity coordinates. The 
results are given by 


■ +(m^u^/HT.)sec^x]secx 

-»■ (2tt)' '^^[1 - 2 ( - q^«pQ/nKTj )^^^exp(q^«PQ/HTj) 

+ erf ( ( - q^flPo/**'*^j ) 3 ’ ^ 

+ (- 2q^<pQ/HTj)^''^(mgU^/KT^)’''^secX 

-»-3[l - (q^qpo/xTj iJCm^u^nTj) sec^x(exp(q^Q)Q/KTj) 

(6a) 


and 
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VNj">e> - 

(2tt)' '^^[l - 2 ( - qjfPo/nKTj ( ’^*exp (q^p/HTj) 

+ erfU-qgiPo/HTj)'''^)]-' 

X [{2/3) (- 2qgq)o/«T.)^''^(mguV’>Tj) *'^secx 

■*■ ( ^ +2] exp(qgq)g/KTj 


- (digU^/HTj ) ^^^secTC. (6b) 

Typical values of normalized exospheric heat fluxes are plotted 
in Figure 2a, which includes the indi'idual components Qp ^ 
and the consolidated elements Qn and TIQp \ 

and the total QsQj|+Q, for conditions ( » 0.06Mj^, T 2 * 6T ^ ) 
described as typical by Feldman et al . (I97^a). Figure 2b pro- 
vides normalized plots of the total Q for other values of T 2 /T^ 
and ^2/^p* 

A numerical example would help to illustrate the use of 
Figures 1 and 2. For parameters (dundhausen , 1970) characteris- 
tic of the quiet solar Wind (u = 320 km/sec, = 1,2 x 10^ ®K) 
one obtains secx = 1,61 and thus (2Trm^u /«T^) ' secx = 0,957. 

It follows from Figure la that - 1.^84, which is 

to say that cpQa-15.4 V for N 2 = 0,06Mp and T 2 /T^ * 6, Locating 
this value of ( " PgCpQ / *^T ^ ) In Figure 2. one obtains (for ex- 
role) the result that ( m^/nT, ^ (0/N^m_ ) = 2.257, which is 

6 1 P ^ 


- il- 








- 2 2 - 3 

to say that Q =» 4,04 x 10 erg/cm sec for « 8 cm the 

proton density recommended by Hundhausen (1970). This heat 

flux should be compared with the corrected value 
-2 2 

(Q ^ 1-3x10 erg/cm sec) obtained observatlonally by Montgomery 
et al . (1968). However, using the N 2 /Np * 0 contour in Figure 1, 
one would have obtained “ 1.094 or cpQ*te11,3 V (cf . 

Schulz and Evlatar, 1972). Using (6) to evaluate Q In this 
case (which Is included In Figure 2b) one obtains 
( m^ /»tT I )^^^( Q/Npme) = 0,271 or Qa«0.72x 10 ^ erg/cm^sec for 
Np » 8 cm ^ and = 1,2xlO^®K. The hot component (j =2) 
thus contributes most of the electron heat flux at r = rQ in 
this fully exospheric model. 
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3. EXOSPHERIC TREATMENT OF HOT COMPONENT ONLY 


The ideal electron exosphere described above is unstable 
to a variety of plasma wave modes (Forslund, 1970; Schulz ana 
Eviatar, 1972; Perkins, 1973). Specific predictions have been 
muddled because of an unwillingness of authors ( e.a . . Schulz 
and Eviatar, 1972) to do high-beta calculations for all of the 
relevant wave modes, but current opinion ( e.g. . Hollweg, 197^) 
favors the obliquely propagating electromagnetic proton-cyclotron 
mode as the dominant instability, operative wherever the peaks 
(maxima) of the proton and electron distribution functions are 
separated by ~ c^ (the Alfven speed) or more, in velocity space. 
For the quiet solar-wind parameters proposed by Hundhausen (1970) 
one obtains u sec x as 51 5 km/sec and c^w38.7 km/sec, 

Hollweg (197^) conjectures that the resulting instability 
reduces the veloc i ty- space separation between the electron and 
proton peaks to a value v^wO.230p “^0 ^ c^, where gj » BrrNj kT^ /B'" 

for protons (j = p) and electrons (j = e), essentially by ad- 
ding usecX-v^ to the parallel velocity of each electron In 
the distribution. Various arguments have been given for and 
against this conjecture. However, the observations of Feldman 
et al . (1974a) seem to bear it out, as applied to the major 
(cooler) electron component. In other words, the majority of 
solar-wind electrons belong to a distribution having ih» form 

r)«{mg/2TTHT ^ ' - m^vJ/axT 

2 

x exp [- ( m^/2vT V|, - u sec X ) ], (7) 
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(with -n 1 } as ar) empirical fit to the observational data. 

This distribution is not exospheric, since it is not an even 
function of v„ for electrons satisfying (1), nor does it vanish 
exactly for electrons ^ailing to satisfy (1) at negative v^. 

However, it may still make sense to treat the minor (hot) 
electron component exospherical ly , as a first approximation. 

On doing so, one obtains 

( 2 nm^uVHT 2 )^^^sec x + 

« 2[1 - (q^cppM^ )]exp(q^(PQ/HT2 ) 
x[l+erf((-q^<Pj,AT2)^^^) 

instead of (4). Comparison of (4) and (8) reveals that qe«Po/*‘^2' 
is the same function of ( 2 TTm^u^/xT 2 )'^secx + ( 2 nmgri^c^/HT 2 )'^^(N^/N 2 ) 
in (8) as (* qgqpo/’^T^ 1 ) 'S of ( 2TrmgU^/KT ^'^^sec x for N 2 /N « 0 

r' 

in (4), Thus, the potential cpQ given by (8) can be read from 
the contour = 0 in Figure I if the ordinate and abscissa 

are appropriately relabeled. Continuing the numerical example 
of the previous section ( N 2 ■ 0,06Np, T 2 « 6T^ « 7.2 x 10^ °K, 
u ■« 320 km/sec, c^ ■ 38.7 km/sec), one obtains a value of 0.9037 
on the abscissa for n ■ 1 . The corresponding value of {-qgCpg/HT 2 ) 
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on the ordinate is 1.169, which is to say that qpQw72.6 V. 

The electrostatic potential aPq is significantly larger in this 
case than for the purely exospheric treatment of both electron 
components. 

( 2 ) ( 2 ) 

The exospheric heat fluxes Qi) and remain as given 

algebraically by (6), but the relation between ( ’ 9 e«pQAT 2 ) 

2 1/2 

and (mgU/’'T 2 ) sec x therein must be obtained from (8) rather 
than from (4). Moreover, the heat fluxes carried by f vy , Vj^; rg) 
are now given by 

= • (l/2)[3 + (mgT,^cJ/HT,)](tn^,n^c|/HT,)’'^ (9a) 

and 

VN,m^) - - ( (9b) 

3/2 

rather than by (6). The various contributions to (m^/HT^'^' 

(Q/N m ) have been plotted (using t\ w > ) In Figure 3a for com- 
P ® 

parison with Figure 2a. Inserting opQ = 72.6 V and = 8 cm'^, 
one obtains ( mg/*<T Q/Npm^ ) = 1.519 or Qw 2.72 x lO’^erg/cm^sec. 
This heat flux is considerably smaller than that obtained above 
In the purely exospheric two- component model. Moreover, the 
total electron heat flux is essentially that of the hot com- 
ponent, as Feldman et al . (1974a) have found by analyzing the 
observed distribution function f(vn,Vj^;r) at r - 1 All. 
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(mgi72c^-//cTj)'^2 = 0.064 









The contours plotted In Figure 3b correspond more nearly 
to as would be required (Schulz and Evlatar, 1972) for 

maqnetoson ic Instability at r ^ 1 AU, For the solar-wind pa- 
rameters used above, this choice of r\ corresponds to an abscissa 
of 1.417 (and therefore an ordinate of 0.6122) on the ^ 

contour In Figure la. This value on the ordinate Indicates 
an electrostatic potential cPq^ 38.0 V and (In Figure 3b) a nor- 
malized heat flux ( mg/*<T ^ Q/N^m^ ) = 1.765; for Np * 8 cm 

- 2 2 

the corresponding physical heat flux Is Q«3*16xlO erg/cm sec. 

- 2 2 

This lies between the value for -nttel (Qw2. 72x10 erg/cm sec, 

- 2 2 

as shown above) and the exospheric value (Q«4.04x 10 erg/cm sec). 
Similarly, the electrostatic potential obtained for « 2 (epQ»38.0 
V) lies between that obtained for n w 1 (cpqsw 72,6 V) and that ob- 
tained for the exospheric case (?Pqw 15.4 V). 
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. DISCUSSION 


The foregoing results suggest the subtle manner In which 
plasma Instabilities may reduce the heat flux carried by solar- 
wind electrons belonging to the hot (j » 2 ) component. Although 
(In this model) only the cooler (j - I ) electron population 
interacts with the unstable waves, the consequently required 
Increase In makes the quasi -exospher Ic hot component less 
skewed In velocity space. The result Is a diminished heat flow 
moment, as given by" (6). Thia Interpretation emphasizes the 
macroscopic (non-local) character of the phenomenon. 

In the case of a uniform plasma, it might be argued that 
only a small portion of even the cooler electron velocity dis- 
tribution could participate In the resonant wave-particle inter- 
actions. Thus, It would be difficult to understand how the 
cooler distribution as a whole could acquire a mean velocity 
— u secx-T|c^. in the Interplanetary medium, however, the un- 
perturbed motion of the typical electron is a bounded oscil- 
litlon between a magnetic mirror point (r ■ r^) and an electro- 
static mirror point (r - r 2 >r^). The particle goes through 
resonance twice in traveling from r^ to r^, and receives an 
impulse from the unstable waves on each occasion. Involvement 
of the entire cooler electron distribution becomes plausible 
when viewed In this way. Hoi Iweg (1974) conjectures that the 
distribution thus acquires a mean velocity ^usecx-nc^ In the 
frame of reference that corotates with the Sun. The observations 
of Feldman et al . (1574a) confirm t» is conjecture. 

The present model does not quite account for the observed 
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form (Feldman et al. . 1974a) of the hot-electron distribution 
f 2 (v|,v^;r). This is not a major cause for concern. A model 
that real i stical 1y incorporates collisions and wave-particle 
interactions will presumably yield a more nearly Maxwellian 
f 2 (v r ) as observed, without drastically modifying the electro- 
static potential cpQ, The vigilant critic might ask why the 
base of the exosphere (r =* rg) should happen to coincide with 
the position of the observer (r ~ 1 AU), as is tacitly assumed 
in applying (l)-(9) to local plasma conditions. This question 
serves to emphasize the quasi -exospher ic spirit of the present 
work and similar past works (e.g, , Perkins, 1973; Hollweg, 1974). 
The plasma is being treated locally as i f it were at the base 
of a true exosphere. However, the exosphere itself plays tha 
role of a "straw man", corresponding to the role of Vlasov equi- 
libria if) simpler plasma-kinetic problems. 

In the limiting case Nj/Np * 0, the results summarized 
by (4) and (6) reduce to the expressions obtained by Hollweg 
(1974). The present results for cpQ bear a close resemblance 
to those obtained by Feldman et al. (1974b) from a more nearly 
hydrodynamica’ analysis. The quas i -exospher ic treatment of- 
fered here should be considered as a viable alternative to the 
fluid phi losophy. 

The electrostatic potential cpq - 40-70 V, obtained by in- 
serting local solar-wind parameters in (8), corresponds to only 
a modest increase (<2 km/sec) in the proton bulk velocity u 
between r ~ 1 AU and r ■ », However, the gravitational potential 
of the Sun amounts to only -9.4 eV/nucleon at r - 1 AU. If the 
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electrostatic potential scales as 1/r, at least back to r - 10-20 
(solar radii), then the corresponding radial electric field 
is ample to overcome gravity and accelerate solar-wind protons 
to the observed bulk speed u - 320 km/sec. Such exospheric ac- 
celeration would correspond qualitatively to the observed out- 
ward proton heat- flow moment (Schulz and Eviatar, 1973). Alpha 
particles in the same electric field would experience a smaller 
acceleration of the same sign, since the charge/mass ratio is 
half that of protons. It is difficult to understand^ j n this con- 
text, the occasional observations by Asbridge et al . (1974) 
that alpha particles have slightly I araer bulk velocities than 
protons in the solar wind. 

In summary, quas i -exospher ic models help to clarify cer- 
tain observable properties of the interplanetary plasma and 
serve as a point of departure for the Vlasov analysis of plasma 
turbulence in the solar wind. Such models are strikingly in- 
complete, how^ zr, and leave many important questions unanswered. 
If the method of quasi linear theory could be adapted to the 
geometry of this inhomogeneous and almost col 1 ision less medium, 
the role of wavelike turbulence in the interplanetary plasma 
might be ascertained In quantitative detail. 
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atmospheric pollution, and high- power gas lasers. 
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pheric~optirs. chemical reactions in polluted atmospheres, chemi' al reactions 
of excited species in rock.-t plumes, chemical thermodynamics, plasma and 
laser-induced reactions, laser chemistry, propulsion chemistry, space vacuum 
and radiation effects on msterials, lubrication and surface phenomena, photo- 
sensitive materials and sensors, high precision laser ranging, and the appli- 
cation of physics and chemistrv to problems of law enforcement and bicmedicine. 

Electronics Research Laboratory : Electromagnetic theory, devices, and 
propagation phenomena, including plasma electromagnetics: quantum electronics, 
lasers, and electro-optics: communication sciences, applied electronics, semi- 
conducting, superconducting, and crystal device physics, optical and acoustical 
imaging; atmospheric pollution: millimeter wave and far-infrared technology. 

Materials Sciences Laboratory : Development of new materials; metal 
matrix composites and new forms of carbon; test and evaluation of graphite 
and ceramics in reentry: spacecraft materials and electronic components in 
nuclear weapons environment; application of fracture mechanics to stress cor- 
rosion and fatigue -induced fractures in structural meta's. 

Space Physics Laboratory : Atmospheric and ionospheric physics, radia- 
tion from the atmosphere, density and composition of the atmosphere, aurorae 
and airglow: magnetospheric physics, cosmic rays, generation and propagation 
of plasma waves in the magnetosphere; solar physics, studies of solar magnetic 
fields; space astr '>my, x-ray astronomy; the effects of nuclear explosions, 
magnetic storms, ai.d solar activity on the earth's atmosphere, ionosphere, and 
magnetosphere; the effects of optical, electromagnetic, and particulate radia- 
tions in space on space systems. 
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